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ABSTRACT 


An investigation into the feasibility of using a 
solid fuel, dual combustion mode ramjet in a high Mach 
number tactical missile was conducted to determine the 
regimes where this propulsion plant might be superior to a 
conventional solid fuel ramjet. The results of the 
study show that at Mach 6, 80,000 feet, with a temperature 
rise combustion efficiency of 90%, the DMRJ performance 
exceeded that of the SFRJ by 20% at near stoichiometric 
overall fuel-to-air ratios, neglecting mixing losses in the 
supersonic combustor. It also appears that nay must be 
greater than 70% in order for the DMRJ to outperform the 


conventional SFRJ. 
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I. INTRODUCTION 


The purpose of this study was to investigate the 
feasibility of utilizing a solid fuel, dual combustion mode 
ramjet (DMRJ) as the propulsion plant for a high Mach 
number tactical missile. As such, the primary objective 
was to examine the conditions under which the performance 
of a DMRJ might be superior.to that of a conventional solid 
fuel ramjet (SFRJ) at a Mach number of 6.0 using realistic 
propulsion performance parameters. The secondary objective 
was to experimentally examine the feasibility of sustaining 
combustion in a Mach 2.5 flow. 

There are many potential benefits from the development 
of a hypersonic vehicle. Among them are the nation's 
efforts to construct the aerospace plane for hypersonic 
travel through the upper and lower atmosphere for both 
Civil aviation as well as space exploration. Benefits to 
the military include powerful tactical/strategic missiles 
and interceptor/fighter aircraft. 

Past tactical missile designs have relied almost 
exclusively on solid rocket motors. Their simplicity of 
design, lightweight construction and ease of handling and 
storage makes them ideally suited for the environment 
aboard ships. Rockets have very high thrust-to-weight and 
can deliver reasonably high specific impulse. They are, 
however, limited to short ranges in order to have suffi- 


cient energy during the end-game of the intercept to 
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defeat defensive maneuvers by the target. The Navy's 
current suite of air-to-air missiles include the AIM-7M 
Sparrow, the AIM-9M Sidewinder and the AIM-54C Phoenix. 
All of these weapons have relatively short range with high 
thrust-to-weights, capable of attaining speeds on the order 
of Mach 1.5 to 2.5 above launch Mach. 

With the advent of very long range detection radars, 
satellite targeting, over-the-horizon data link targeting 
and long range cruise missiles, these weapons are fast 
becoming incapable of adequately defending the fleet in a 
tactical scenario. New weapons like the Tomahawk and 
Harpoon cruise missiles utilize turbojet propulsion that 
has allowed greatly extended ranges. However they were 
designed for slow, non-maneuvering surface targets and are 
themselves subsonic. They are, therefore, practically 
useless against a tactical threat. What is needed is a 
long range, standoff weapon capable of reaching out 30 to 
100 nautical miles with a high terminal velocity giving it 
sufficient energy to defeat the threat. 

In order to achieve long ranges and high velocities, 
specifically terminal velocity, the vehicle must be powered 
Chreoughout the full range of the flight. Rockets use boost 
and sustain propellant grains that burn out prior to 
intercept, generally leaving the missile to ultimately 
coast to the target. A pulsed motor is one means of 
extending the range of solid rockets. Another viable 
propulsion device for longer ranges is a high speed, 


airbreathing engine. Below Mach 2.5, the gas turbine 


family of engines is inherently more efficient and will 
outperform all other air breathing power plants. At speeds 
greater than Mach 3.0, however, ram pressure alone is 
sufficient to provide proper compression of the airflow so 
that the compressor actually becomes a liability [Ref. 1: 
p. 38). In addition, materials limit the maximum tempera- 
ture allowed in the turbine section so that as Mach number 
increases, less and less energy can be added to the flow 
before that limit is reached. Thus, thigusieeealls oft 
rapidly on the high speed end for a turbojet. 

Between Mach 3.0 and 5.0, the ramjet becomes more 
efficient than the turbojet. Although ramjets have been 
around for some time, it has been only recently that their 
primary disadvantage has been overcome. Below Mach 0.9, 
ram pressure is insufficient to compress the flow adequate- 
ly for efficient combustion [Ref. 2: p. 140]. Thus, a 
ramjet must be brought up to speed before it can be 
employed. The advent of the integral-rocket-ramjet (IRR) 
has made this possible, by employing a solid propellant 
booster to take the ramjet up to its take-over Mach number. 
Figure 1 presents the three power plants discussed thus 
far. } 

As the Mach number increases towards Mach 5.0, shock 
compression to subsonic velocities results in large losses 
in stagnation pressure. Here also, a temperature limit is 
reached for the combustion process due to dissociation of 
reactants in the flow. In addition, since subsonic 


combustion is most efficient when the incoming Mach number 
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Figure 1 - Propulsion systems for tactical missiles: 
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is on the order of 0.2 to 0.3, the static temperature and 
pressure in the combustor approach material limits. 
Figure 2 delineates the limits of the typical ramjet flight 
envelope. 

In order to achieve hypersonic velocities, scramjet 
technology is the most feasible solution. The decrease in 
seeks losses (reducing the combustor flow only to low 
supersonic velocities) more than offsets the increased 
heat addition losses from adding energy to a supersonic 
flow. The lower static temperature associated with 
Supersonic flow ensures very little dissociation will 
occur and the need for cooling is reduced. This, in 
conjunction with lower pressures, reduces the structural 
limitation on the combustor itself. Thus, greater energy 
addition is possible, resulting in higher exhaust veloci- 
ties and more thrust. 

There are two types of scramjets. The first is a 
conventional ramjet with supersonic flow pneu the 
combustion chamber. This is accomplished by designing an 
inlet without a normal shock so that the flow remains 
supersonic. The second type of scramjet ‘is more appro-- 
priately named a dual-combustion mode ramjet (DMRJ). This 
configuration involves using a conventional ramjet as a gas 
generator that exhausts into a second stream of supersonic 
bypass) air. Mixing and combustion of the fuel-rich 
effluent occurs, enabling further heat release. Figure 3 
illustrates both types of scramjets. The primary problem 


with the combustion process is whether or not sufficient 
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Figure 3 - Supersonic combustion systems for tactical 
missiles. [Adapted from Ref. 2: p. 142] 
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residence time and mixing exist for sustained combustion, 
given the speed of the flow and the kinetic rates of the 
process. 

Figure 4 illustrates the ideal performance of the 
propulsion devices discussed thus far, and the flight 
regimes in which they have advantages over each other. 

There are several advantages associated with the use 
of solid fuels in both rockets and ramjets. Liquid fueled 
ramjets (LFRJ) require a positive fuel control, a means to 
atomize the fuel droplets and a means of stabilizing the 
flame. These requirements make the the LFRJ inherently 
more complex, as well as heavier, than the SFRJ. 

The SFRJ is "“self-throttling". The amount of mass 
flow taken into the engine and its properties determines 
the fuel flow rate in the combustor. The general expres- 


Sion for the regression rate for an SFRJ is: 
(eqn. 1) r= G* pY té4 


where G represents the mass flux entering the combustor, p 
is the chamber pressure, t is the air inlet (static) 
temperature and x, y and 2 are specific to the fuel 
and generally lie between 0.1 and 0.8. The flame stabili- 
zation is accomplished» through the geometry of the combus- 
tor rather than insertion of vee-gutters or other bluff 
body flameholders. Figure 5 illustrates the flow field in 


the flame stabilization region of an SFRJ combustor. 
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Figure 4 - Theoretical performance envelopes. 
[Adapted from Ref. 2: p. 144] 
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mlchough the SFRJ is "self<-throttling", the throttling 
is not optimum. The positive fuel controller in the LFRJ 
Will enable it to achieve higher performance throughout the 
flight envelope than the SFRJ. This situation is changing 
as metallized fuels enable a higher energy density in the 
SFRJ fuel grain (resulting in improved performance) and 
variable air flow geometries are being developed. An SFRJ 
is able to achieve very high performance with a simple 


design. 
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II. NATURE OF THE PROBLEM 


A. ANALYTICAL SOLUTION: SOLID FUEL RAMJET 


For the purposes of this study, analytical models had 


to be assumed for both the SFRJ and the DMRJ. For the 


SFRJ the following conditions were employed: 


(a) 


(b) 


(c) 
(4) 
(e) 
(f) 
(9) 
(h) 
(1) 

(4) 
(x) 
(1) 
(m) 


Design point: Mg = 6.0; h = 80,000 ft; 1 lbm/sec 
flow rate 

Four-shock inlet with ramp angles of 12, 28 and 43 
degrees (from horizontal), or an ideal inlet with 
no shock losses 

Inlet operating 10% supercritical 
Subsonic pressure recovery = 0.95 in the inlet 
Combustor inlet Mach = 0.2 i 

Constant area combustor 

One-dimensional flow in the combustor and nozzle 
HTPB grain as fuel (p = 0.0332 lbm/in?) 

Cp = 1.5 in the SFRJ combustor 

Equilibrium, adiabatic Sonnoeticn 
Nozzle throat choked 

Isentropic exhaust nozzle 

Frozen flow expansion through the exhaust nozzle 


to ambient pressure 


Figure 6 depicts the station numbers in effect for the 


analysis. 
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Figure 6 —- Solid fuel ramjet model. 
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The continuity equation (conservation of mass), in the 


form below, was used to size the cowl and the combustor. 











ont 
P v—( 2 i 
RT\y+1 
(eqn. 2) m= 
A/A 
M 


Oblique and normal shock relationships were used _ to 
determine the pressure loss up to station 2. With Mach 
number and flow rate specified and stagnation temperature 
known, equation 2 yields the area, Ap2. } 

The only known parameters after the combustor were the 
throat Mach number and the exit pressure. To proceed 
across the combustor, conservation of energy and momentum 
are applied. For the energy equation, the Naval Weapons 
Center (NWC) Propellant Evaluation Program, PEPCODE, was 
used for equilibrium, adiabatic combustion to construct 
tables of combustion gas properties and theoretical 
stagnation temperatures at station 4 as a function of the 
stagnation properties entering the combustor and the 


fuel-to-~air ratio (f) [Ref. 3]: 
(eqn. 3) PEPCODE = fn(f,T,P) => Tin, Cpe Ry 


Three tables are included in Appendix A. In order to 
determine the properties at station 4, knowing the proper- 
ties from PEPCODE, the following equations had to be 


simultaneously satisfied to determine M,: 
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(eqn. 4) 
(eqn. 5) 
where 
(eqn. 5a) 
(eqn. 6) 
(eqn. 7) 
(eqn. 8) 
(eqn. 8a) 
(eqn. 9) 
RAMJET, 


solve this system of equations, given the output of PEPCODE 
and an initial guess for M, as the input data. 
iteration was checked with the continuity equation until a 
specified tolerance was met. 
was 0.0001 in mass flow rate. 


and specific fuel consumption (SFC) were then determined 
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an interactive computer code, 
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Was MW Gl etcet tO 


In this case, that tolerance 


The exhaust velocity, thrust 


y-1 





P,. 
(eqn. 10) uv 2Cg,s7,|1- (=) Y 
e 
mute—om~-miu 
(eqn. 11) F= Ch ee 
ee 
i 
(eqn. 12) SFC = — 
F 
A sample output is included in Appendix B, along with 


the progran. 


The grain length is determined by specifying 6 for each 


case. Equation 1 determined the regression rate with, 


« = 7.1789x107° 
xX = 0.53 
y = 0.33 
z= 0.71 
for HTPB. With fuel flow rate specified, the grain 
length, Lg, is determined from, 
(eqn. 13) m,=pmD Lor 


B. ANALYTICAL SOLUTION: DUAL COMBUSTION MODE RAMJET 

As with the SFRJ, a model had to be created to deter- 
mine all the properties through the DMRJ. In their 1958 
paper on supersonic combustion, Weber and MacKay [Ref. 5] 
studied a one-dimensional model of a pure scramjet with a 
constant area combustor. Their model used liquid hydrogen 
as the fuel, for which a fixed combustion efficiency of 95% 


was assumed. They accounted for neither momentum losses 
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due to fuel injection in a high speed flow and shock losses 
in the combustor, nor were wall friction losses applied in 
the combustor. Dissociation was addressed briefly. For 
their study, the scramjet outperformed the conventional 
ramjet over the entire range of combustion temperatures 
examined (taken at an altitude of 120,000 ft. Mach 6) as 
long as there was no recombination in the exhaust nozzle 
expansion process (i.e., frozen flow) [Ref. 5: p. 106,110]. 
This stresses the impact on performance that dissociation 
has due to very high static temperature in the ramjet 
combustor. For an isentropic inlet, their scramjet model 
was superior to the conventional ramjet above Mach 7. With 
a wedge-type inlet, that advantage shifted to Mach 5. 
mee. Sp. 102] 

In this study, an attempt was made to take advantage 
of known realistic performance parameters for momentum loss 
in the SFRJ combustor and known subsonic combustion 
efficiency behavior of the solid fuel, HTPB. A one-dimen- 
Sional, constant area combustor was also employed. Equili- 
brium adiabatic combustion through PEPCODE accounted for 
some of the effects of dissociation, though not completely, 
and yielded variable properties. Rather than assume a 
value of combustion efficiency, n,7 was used as a parameter 
in the performance study. 

Figure 7 depicts a schematic representation of the DMRJ 
used in this study with applicable station numbers. 
Upstream of the "B" path of the mixer was an SFRJ gas 


generator identical to Figure 6. For use with the DMRJ, 
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Figure 7 - Dual-mode ramjet mixer model. 
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however, station 4 has been renumbered as’ station 3. 
Two separate inlets were assumed; one for the SFRJ gas 
generator and one for bypass air. All the assumptions from 
the previous section for the SFRJ were applicable to the 
gas generator except in this case, the total design mass 
flow rate was 1 lbm/sec. For the bypass portion, the 
following conditions were employed: 

(a) Three-oblique-shock inlet with the identical 
pressure recovery as the SFRJ inlet (prior to the 
normal shock), and an ideal inlet with no losses 

(6) One-dimensional flow 

(c) Mixer inlet Mach = 2.5 

(ad) Constant area supersonic mixer 

(e) No mixing losses 

(f) Static pressure at the mixer inlet determines 

exhaust pressure for the gas generator exhaust 
(i.e. SFRJ nozzle sized to match pressures at 
station 3) 

(g) Equilibrium adiabatic combustion 

(h) Frozen flow expansion to atmospheric conditions 

(1) Isentropic nozzle 

The method for solving the SFRJ portion of the DMRJ 
was identical to the previous section, except that the 
exhaust pressure was now equal to p3a, vice atmospheric. 
RAMJET was again used to determine all the properties at 
station 3B. As before, oblique shock relationships were 
used to determine the properties at station 3A. Equation 2 


was then used to determine the areas at stations 3A and 3B. 
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The mixing process presented the most difficulty in 
the analysis. No information was available in the open 
literature to make a reasonable assumption for the mixing 
loss. The approach taken was to assume a value of thermal 
efficiency (nar), representing losses in the heat release 
process only, to vary the equivalence and bypass ratios 
and to not account for the mixing loss at all. From this,’ 
a "baseline" performance could be identified to which 
mixing losses could be applied as they become known. 

PEPCODE was again employed to determine the properties 
at the exit of the mixing chamber, given the stagnation 
temperature of the Mach 6.0 flow, an estimate of the 
chamber stagnation pressure (to be verified) and the 
fuel-to-air ratio (f). In this analysis, a mass averaged 
temperature was determined for station 3, assuming complete 


mixing and no combustion, defined as: 


(eqn. 14) ie = MaaTpaaT aa * apCpaa7 an 
(m3, + m,)C, 
Proceeding as before, with eta specified arbitrarily, all 
the properties known at station 3, and station 4 gas 
properties from PEPCODE, all the flow properties at station 
4 were determined by solving equations 2,4,6,7 and 9 
simultaneously. 

MIXER, another interactive code, was written to solve 
this simultaneous set. Input variables included all the 
properties at stations 3A and 3B, the PEPCODE output for 


station 4 and an estimate of the Mach number at station 4. 
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Rather than using the continuity equation as the check, 
this program determines the velocity of the flow at 
station 4 by two means: (1) by equations 4,6 and 7; 
(2) by conservation of momentum, equation 8. Repeated 
estimates of the Mach number are made until the difference 
between the two velocities were on the order of 0.1% 
This program and a sample output are included as Appen- 
dax Cc. 

Once the properties at station 4 were determined, the 
exhaust velocity, thrust and specific fuel consumption 
were again determined using equations 10 - 12. This 
calculation was performed for each bypass ratio considered 


and each value of thermal efficiency. 
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IIIT. ANALYTICAL RESULTS 


A. SOLID FUEL RAMJET 

The performance of the SFRJ is shown in Figure 8. The 
cases considered were: 1) design Mach number of 6.0, 
4-shock inlet; 2) design Mach number of 6.0, isentropic 
inlet; and 3) design Mach number of 4.0, 4-shock inlet. 
All cases were done for an altitude of 80,000 feet. Tables 
of properties for these engines are included in Appendix D. 

The shape of the curves is highly dependent on the 
adiabatic combustion efficiency, n,;. Figure 9 represents 


a plot of equation 5, which is the efficiency behavior of 


HTPB. The minimum at © = 1.0 explains the pronounced dip 
in the performance curve. This dip is present with or 
without losses, regardless of design Mach number. The 


lower design Mach number, however, is more sensitive to 
changes in nar. For example, for an overall equivalence 
ratio of 1.3 (nap = 0.75), a 1% change in efficiency at 
Mach 4 will cause a 37% greater temperature rise than for 


Mach 6:3 


a ae _1640_ | 


I — — — 
A hie 3011 |4 2200" | 


This temperature rise is directly related to the exhaust 
velocity, therefore the thrust. Thus, trends in the 


performance curve are accentuated at the lower Mach number. 
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The performance curves indicate that specific thrust 
rises with @ despite the behavior of Ty+,, which is shown 
mee rigure 10, The reason for this is two-fold. The 
primary effect is that of decreasing molecular weight. 
This is reflected in higher values of the gas constant, R, 
for higher values of 6 (refer to tables in Appendix A). 
Equation 10 shows that the dependence of exhaust velocity, 
and therefore, thrust, is /R. In addition, increasing 
® implies an increase in mg. Thus, both of the terms in 
the jet thrust (m.u,) increase while the ram drag term 
(MoUg) remains constant. 

The secondary effect is that of the increasing combus- 
tion efficiency of HTPB in the fuel-rich regime (refer to 
Figure 9). The combination of these two effects appear to 
offset the sharp decrease in theoretical adiabatic combus- 
tion temperature, enabling the thrust to increase. 

Inlet losses also have a very iaeae effect on the 
performance of the SFRJ. Figure 8 shows that for the same 
design Mach number, inlet losses reduce the specific thrust 
of the SFRJ by approximately 20% throughout the range of 
equivalence ratios. This is due solely to the stagnation 
pressure drop across the inlet system. 

The Mach 4 curve lies above the Mach 6 curve, as 
expected for a ramjet. The ram drag is the primary 


reason for this, as illustrated below for 6 = 1.8. 
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Figure 10 - Equilibrium adiabatic combustion tempera-~ 
ture for HTPB for isentropic inlet (640 psi) and for 
4-shock inlet (120 psi). 
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Although there is a large increase in exhaust velocity 
in going to Mach 6, there is an even Breates increase in 
the freestream velocity, which leads to higher thrust at 
Mach a 

A second contributor to the higher performance at the 
lower Mach number is the decrease in shock losses due to 
the weaker shocks at Mach 4. 

Dissociation effects are accounted for by PEPCODE. As 
expected for the SFRJ, an increase in design Mach number 
diminishes the effect of adding fuel to the flow. For 
example, for My = 2.8, the difference in theoretical 
stagnation temperature rise across the combustor going 
from @ = 0.5 to 6 = 1.0 (i.e., doubling the fuel flow) is 
approximately 1300 R. For Mg = 6:0, however, the same 
increase in fuel flow leads to a difference in theoretical 
stagnation temperature of approximately 580 R. This is the 
effect of Bee eeeon on energy addition potential. 

The unique combustion efficiency behavior for the SFRJ 
results in a predicted F/mg versus 6 behavior that is 
somewhat different than the behavior that is usually 
reported for ramjets with fixed combustion efficiency. 
Specific thrust will maximize at very rich equivalence 


ratios. 


B. DUAL-MODE RAMJET 
The results from the analytical model support the 
previous theoretical works of both Billig [Ref. 2] and 


Weber & MacKay [Ref. 5]. The performance of the DMRJ 
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exceeds the performance of the conventional SFRJ, within 
the constraints of the conditions applied, by approximately 
20% in the range of anticipated equivalence ratios (0.8 - 
irate) Each parameter was examined for its impact on the 
calculations and results. However, because no attempt was 
made to account for mixing/shock/wall friction losses in 
the supersonic combustor, these results identify a minimum 
required performance, or a "baseline", to which these 
losses must be applied. 

The parameters examined included the combustion 
efficiency (n,;), the bypass ratiol, the design Mach number 
and the overall equivalence ratio. The effects of using 
an isentropic inlet were looked at briefly, as well as the 
impact of dissociation in the supersonic combustor. The 
following bypass ratios were used in the analysis: 50/50 
(1.e., 50% primary air, 50% bypass air), 60/40, 75/25 and 
80/20. Tables of properties through the engines are 
included in Appendix E and will be referred to below. 

It must be emphasized that each point represents a 
different design condition and, therefore, a different 
engine. Each engine was re-sized as a parameter was 
changed and then flown at its design point. (The design 
point was Mach 6, 80,000 feet unless otherwise specified.) 

1. Performance Comparison: nar 

As with the SFRJ, the performance of the DMRJ is 


highly sensitive to the combustion efficiency specified. 


. ES Sie ae is defined as the ratio of mass flow 
through the primary gas generator to the mass flow through 
the bypass duct. 
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The gas generator has the same efficiency characteristics 
as the SFRJ (refer to Figure 9). The efficiency in the 
supersonic combustor was varied between 0.25 and 0.9. 

Figure 11 is a plot of specific thrust versus the 
overall equivalence ratio for the conventional SFRJ, the 
DMRJ with an efficiency of 90% for all bypass ratios and 
the DMRJ for a combustion efficiency of 25% for two 
bracketing bypass ratios. 

To understand the effect that efficiencies have on 
the DMRJ, it was necessary to compare an isentropic SFRJ to 
an isentropic DMRJ with nar = 1.0 for all combustion pro- 
cesses. In this case, the engines were isentropic in that 
there were no inlet losses present. For the table below, 
Mp = 6.0, h = 80,000 ft. and Smop = 1.0 (1.€e., same total 
fuel flow). 


SFRJI 80/20 50/50 
sorry —- aod ee D 2.00 
T3 alae 5600 4848 
T4 5565 5565 5565 
Pi 593 445 291 
M4 0.93 2.10 1.48 
Lede | 45 193 347 
T3ave —— 5230 4166 
AT =—— 53.5% 1399 
Ty 5490 3544 4355 
F/Mo 1O37..6 101.4 97.4 


What this information shows is that at Mach 6, 
80,000 ft., the SFRJ has superior performance to the DMRJ. 
The reason behind this is the stagnation pressure drop due 
just to heat addition. For the SFRJ with subsonic heat 


addition, the drop is very small. It, therefore has the 
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Figure 11 - Performance comparison for SFRJ, DMRJ with 
nat = 0.9 for four bypass ratios and DMRJ with nar = 
0.25 for two bracketing bypass ratios. (Myo = 6.0, 


80,000 feet) 
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highest stagnation pressure. For the 50/50 DMRJ, the 
cooling effect of the bypass air is very significant 
(T3ave = 4166 R). With considerable unburned fuel and a 
low inlet temperature, much of the heat addition occurs 
in the supersonic flow. Thus, it has the largest pressure 
axep. The 80/20 DMRJ lies between these two extremes. 
This result supports oe finding of Weber & MacKay that for 
the same heat addition (both types of engines have the same 
T4) the engine with the highest P, will have the highest 
exhaust velocity, which means highest thrust [Ref. 5: 
erro). 

nat Will affect the stagnation temperature at the 
exit of the combustor, supersonic or subsonic. When 
NxaT = 0-9 in the supersonic flow, the DMRJ exhibits a 
Signifigant performance advantage at lower values of 
owerall equivalence ratio. However, on the fuel-rich end, 
the DMRJ levels off while the SFRJ continues to increase, 
eventually causing a cross-over in performance. This 
behavior is explained in the table below for a DMRJ with a 
bypass ratio of 80/20: 


: DMRJ SFRIJ DMRJ SFRI 
pre. TING NEC Ta cael | Oe em a. OO 
on 0.88 0.90 | 1.76 1.80 
T3Bth 5391 5292 | 4581 5031 
T3p 4693 4626 | 4612 4850 
Ty, 5195 4626 | 5008 4850 
F/Mo TalleaS 60.6 | 91.1 91.5 
Nati Ro 0.68 0.68 | toe 0.90 
nap ao 0.90 ~- | 0.90 =< 


The first two columns represent a comparison at the 


fuel-lean end ( mop = 0.9) where the DMRJ has a distinct 
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performance advantage. The third and fourth columns 
compare the SFRJ and DMRJ at the fuel-rich end (mop = 1.8) 
where the performance crosses over. © 

For the first case, the SFRJ efficiencies are 


identical so that the temperature at the exit of the 


combustors are nearly the same. The only difference was 
that @ = 1.1 lies closer to the peak of the temperature 
profile, shown in Figure 10, than does $ = 0.9, leading to 


a higher theoretical stagnation temperature, hence a higher 
actual stagnation temperature. T, is then significantly 
higher due to very efficient burning of the excess fuel in 
the second combustion process. This is the driver behind 
the performance advantage for the DMRJ at the lean end of 
overall equivalence ratio. 

For the second case, the performance is very nearly 
equal, with the SFRJ being slightly superior. In this 
case, the DMRJ primary combustor is very fuel rich at a $ 
of 2.2. This point is well past the peak in theoretical 
stagnation temperature for HTPB. The DMRIJ actually 
exhibits a drop in stagnation temperature at the primary 
burner exit from an increase in the fuel flow while the 
SFRJ exit temperature increases. The efficiency behavior 
of HTPB causes the actual exit temperature to ‘increase 
slightly, but it is still less than the exit temperature 
for the SFRJ. The crossover in performance for the SFRJ 
comes from a combination of higher stagnation pressure and 
smaller molecular weight (larger R) that just offsets the 


difference in Ty. 
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Figure 12 is a plot of DMRI specific thrust versus 
overall © for two bypass ratios with n,; as a parameter. 
This plot shows the increment by which the performance 
varies with efficiency. It appears that for any realistic 
advantage from the dual-mode to exist, combustion efficien- 
cies on the order of 70% should be attained (refer to 
Ruoures 11 & 12). In the general region of interest 
(l.e., &® < 1.2), 50% efficiency leads to near equal 
performance to the SFRJ accounting for heat addition losses 
but no mixing or shock losses. 

2. Performance Comparison: Bypass Ratio 

Both Figure 11 and 12 demonstrate the dramatic 
effect of decreasing bypass ratio as efficiency decreases. 
For very high supersonic combustion efficiency, the 
performance lines for each bypass ratio very nearly 
overlay. However, as that efficiency is decreased, the 
performance lines begin to diverge. For the same value of 
overall ¢, there is a large difference in performance 
between a bypass ratio of 80/20 and one of 50/50. In the 
table below, a closer examination is made for Mp = 6.0, 


80,000 ft. with inlet losses and with nar = 0.5: 


50/50 80/20 
drop 0.90 — 0.88 
Pope 1.80 Leek @ 
nat IRI 0.90 0.68 
M¢ 0.066 0.065 
T3th 5031 5391 
T3 4850 4693 
Taave 4193 4463 
T4 4743 4871 
F/mo 62.6 64.5 
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Figure 12 - Effect of nay on two bypass ratios. 
(Myo = 6.0, 80,000 ft., with inlet losses) 
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For very nearly the same fuel flow, there is a 
difference in performance. Even though the 50/50 DMRJ has 
a higher stagnation temperature exhausting into the mixer 
(T3), due to the higher nar, T3ayve is significantly lower 
due to the larger amount of cooler bypass air (Tp = 3198). 
This causes Ty, for the 50/50 DMRJ to be lower than the 
80/20 DMRJ, leading to a higher specific thrust for the 
larger bypass ratio for nearly the same amount of fuel 
addition. 

Another effect of bypass ratio is that to obtain 
high values of overall ¢, larger bypass ratios must be 
used. For example, to obtain an overall equivalence ratio 
of 1.5 would require égppyz = 1.875 for the 80/20 DMRJ. The 
50/50 DMRJ however, would require a primary combustor 
equivalence ratio of 3.0, for which the efficiency data for 
HTPB is not applicable. Thus, for any specified overall 
equivalence ratio, the 50/50 DMRJ will be much more 
fuel-rich than the 80/20 DMRJ. 

One final observation is that although the combus- 
tion efficiency in the primary combustor is higher in terms 
of temperature rise for the fuel-rich engines (lower 
bypass ratios), there will be more unburned fuel present in. 
the supersonic combustor: For a high n,;7 in the supersonic 
combustor, it doesn't appear to matter, and the performance 
of the two cases should be similar. Figure 12 verifies 
this behavior. For the lower n,, in the supersonic combus- 


tor, that excess fuel is burned less efficiently so that a 
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performance difference should exist. This behavior is 
also depicted in Figure 12. 
3. Performance Comparison: Mo 
Figure 13 shows how decreasing the design Mach 
number to Mach 4 (80,000 ft.) affects the performance of 
the SFRJ and the DMRJ. The bypass ratio of the DMRJ 
is 80/20 and the combustion efficiency is 90%. Inlet 
losses are included. As in the case of the SFRJ, going to 
a lower Mach number improves the performance of the DMRJ 
for the same reasons: lower shock losses and reduced ram 
drag. Figure 14 includes an 80/20 DMRJ with nar = 0.6 to 
show the effect of changing the efficiency and the Mach 
number. 
4. Performance Comparison: Inlet losses 
The effect of using an isentropic inlet on the DMRJ 
with nar = 0.9 is depicted in Figure 15. Just as the 
performance for the SFRJ increases, the DMRJ realizes a 
gain in specific thrust of apoE eae 25% at Mach 6, 
80 ,000 Ser. The reason for the improvement, as before, is 
the absence of stagnation pressure losses due to inlet 
shocks. The SFRJ is included in Figure 15 to show erat the 
isentropic DMRJ gives the best performance. 
5. Performance Comparison: Dissociation 
One of the reasons for going to supersonic combus- 
tion over subsonic is dissociation of the reactant mole- 
cules due to the extremely high static temperatures that 
are encountered when a high speed flow is shocked down to 


low subsonic Mach numbers (0.15 to 0.3). Dissociation 
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MACH NUMBER EFFECT 


85.0 100.0 115.0 


70.0 


LEGEND 


SPECIFIC THRUST (LBF—SEC/LBM) 


95.0 





_ 40.0 


0406 08 10 12 14 416 #18 42.0 2.2 2.4 
PHI 


, 
Figure 13 - Effect of design Mach number on performance 
for SFRJ and 80/20 DMRJ with nar = 0.90. 
(#=cnock inlet, 80,000 ft.) 
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MACH NUMBER EFFECT 


85.0 100.0 115.0 


THRUST (LBF) 


70.0 


LEGEND 


52.0 





40.0 


04 06 08 10 142 14 #16 #18 2.0 2.2 2.4 


5 
Figure 14 - Effect of design Mach number for. SFRJ and 
80/20 DMRJ for two values of nar, 0.60 and 0.90. 
(4-shock inlet, 80,000 ft.) 
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EFFECT OF INLET LOSSES 
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Figure 15 - Effect of using an isentropic inlet as 
compared to a wedge-type inlet. 
(Mp = 6.0, 80,000 ft., nay = 0.9 for DMRJ) 
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prevents complete combustion of the reactants by tying up 
chemical energy in intermediate species. Rather than a 
complete combustion process, where all of the fuel and 
oxidizer react to form CO5 and H:0, species such as OH, 
CO and H are present as well. Their presence reduces the 
potential temperature rise obtained by burning a specified 
amount of fuel and if recombination of the species is 
prevented in the exhaust nozzle, reducing the exhaust 
velocity. This can occur when kinetic rates and residence 
times are such that equilibrium cannot be attained. The 
limiting case is termed "frozen flow". When cinecie rates 
are fast enough compared to residence time in the mere 
or, "shifting equilibrium" exists, where recombination of 
the dissociated species allows more thermal energy to be 
released. 

In using PEPCODE to determine the properties after 
equilibium adiabatic combustion, two approaches can be 
adopted. In this study, stagnation temperature was used 
aS an input variable. This means that the gas properties 
are based on stagnation conditions, i.e., Mg = 0. For the 
SFRI, static and stagnation temperatures are very nearly 
the same so that it is not a problem. For the DMRJ, 
however, this is not the case. My, is supersonic, so that a 
large difference exists between static and stagnation 
properties. A more accurate method would be to use static 
temperature as the input variable in such a way that 
PEPCODE computes properties, based on shifting equili- 


brium, for M, ~ 0. This approach, while more accurate, 
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requires significantly more computation. For each data 
point, PEPCODE must be run twice, as well as RAMJET, to 
solve for the properties at the exit of the combustor. 

This effect was examined and the results are 
presented in Figure 16 for an 80/20 DMRJ, with inlet 
losses, Mach 6, 80,000 ft. and nj, = 0.9. For @< 1, using 
properties based on static conditions has no significant 
effect on the performance of the DMRJ. For 62 1.25, l.e., 
fuel rich, properties based on static conditions would 
result in approximately 5% greater specific thrust. Since 
the anticipated range of equivalence ratios is from 0.8 to 
about 1.1, this effect was neglected and the properties 


were computed using stagnation conditions. 
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0.50 0.75 1.00 1.28 1.50 1.75 2.00 
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Figure 16° - Performance comparison between using 
properties computed by PEPCODE based on static condi- 
tions to properties based on stagnation conditions. 

(Myo = 6.0, 80,000 ft., wedge-type inlet) 
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IV. EXPERIMENTAL PROCEDURE 


The second objective of this study was to assess the 
feasibility of sustaining combustion in a supersonic flow, 
uSing a conventional SFRJ as a gas generator for a superso- 
nic combustion chamber. A schematic of the apparatus is 
depicted in Figure 17. The apparatus was designed and 
constructed, however time constraints prohibited the tests 
from being conducted. 

Two supplies of vitiated air were available, each 
producing a mass flow rate of 0.5 lbm/sec. The design 


parameters are listed below: 


(a) T = 1250 R for both the primary air and bypass air. 

(6b) P = 100 psi for both streams entering the mixer. 

(c) Mach number of 2.5 for both streams coming into the 
Supersonic combustion chamber. 

(4) @= 1.2 in the SFRJ (Lg = 10.59 in.) 


(e) Fuel : HTPB. 


Because the kinetics of the mixing process was largely 
unknown, the streams were joined at an angle of 30 to 
force the mixing. The static pressures of the two streams 
entering the mixer were designed to be nearly equal in an 
effort to minimize losses. Shock losses were expected due 


to the large angle between flows, however these were 
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accepted in order to enhance the chances of mixing and 
combustion. The primary objective, initially, was just 
to obtain combustion, not necessarily efficient combustion. 
The mixer. length was determined from residence time 
considerations. Using an average Mach of 2.0, a mixer 
length of 6 feet should give a residence time of approxi- 
mately 1.3 msec. Temperature and pressure readings would 
be made just entering the mixer and at the exit in order to 
ascertain whether or not sustained combustion occurred. If 
the engine successfully ignited and sustained combustion, 
the mixer would be shortened until combustion eae not be 
attained. This would represent a measure of the kinetic 


rates of the reaction. 


50 





(a) («) 9{[ZZOU BULBIIATP-sUTBAIAUO) 


k—___—_——. 1 9 _______-5 


Lal <a Ss ei 
G) () Das /uqr $*0 


—\/, 


9[ZZOU BUTBIIATP-3ULBIOAUO) 


Figure 17 - Dual-mode combustor apparatus 
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V. CONCLUSIONS AND RECOMMENDATIONS 


From an analytical standpoint, this study showed that a 


solid fuel, dual combustion mode ramjet exhibits areas of 


potential superior performance over a conventional solid 


fuel ramjet for use in a Mach 6 tactical missile. Based on 


the specified flight conditions and performance parameters, 


not accounting for any losses due to mixing in the super- 


sonic combustor, the DMRJ outperforms the SFRJ at near 


stoichiometric equivalence ratios for ma, > 0.5. Addition- 


al results are summarized below: 


(1) 


(2) 


(3) 


(4) 


It is estimated that with mixing losses, combustion 
efficiency in the supersonic combustor must be at 
least 70% or greater to achieve superior 
performance with the DMRJ. 





For supersonic combustion efficiencies less than 
approximately 60%, performance of the DMRJ is very 
sensitive to the amount of airflow in the bypass 


stream. Greater bypass flow led to higher perfor- 


mance for the same overall 96. 


The effect of decreasing inlet losses and lower 
Mach numbers were as expected; both leading to 
better performance. 


Dissociation effects proved to be relatively minor 
( = 1%) in the range of @ near stoichiometric 
( is ~@C. I ® = 0 e 8 = alt ° 1) ° , 


In order to refine the range of required combustion 


efficiencies, follow-on research should focus on the 


actual 


mixing processes, the combustion kinetics and 


supersonic combustion processes. Experimental testing with 


52 


the completed DMRJ apparatus will provide insight into 
some of these areas, aS well as some validation of the 
analytical model. Further studies should be conducted 
using solid fuels with different combustion efficiency 
characteristics than HTPB in order to more closely examine 


fuel effects on the system performance. 
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TABLE 1 
PEPCODE DATA FOR EQUILIBRIUM ADIABATIC COMBUSTION OF 
b if Y R Cy 

0.5 .0367 1.266 53.256 ~326 

0.6 -0440 1.263 53.426 ~330 

Oui7 -0514 1.260 53.689 ~334 

O88 .0587 1.259 54.013 .338 

0.85 .0624 1.258 54.199 ~339 
9 -0661 1.258 54.415 341 
95 -0697 1.258 54.631 ~343 

tO .0734 1.258 54.848 oSa4 

1.05 .0771 1.258 55.095 ~346 

De Ji -0807 1.258 55.342 Noa 

1.15 .0844 1.258 55.605 348 

12 -0881 1.258 55.883 ~350 

1.25 .0918 1.259 56.176 ~351 

3 .0954 1.259 56.470 »352 

19.4 -1028 1.261 £57.119 ~355 

m5 POI 1.263 °° 57.814 .357 

1.6 74a 1.26500858 .5410 »359 

1 .1248 1.268 59.313 361 

Ne 8 1320" 1.275mrCO. YOR ~363 

1.9 ,13905 e1.274—@ 60.878 364 

2.0 -1468 1.277 #£,61.676 . 366 

2 Sai . 280" 62. 450 .367 

22 pres 1.283 63.268 ~ 369 

2053 .1688 1.286 64.056 a0 

Remarks: 


APPENDIX A = PEPCODE Data 


4th te 

4746 1.297 
4938 1.292 
5088 1.289 
5203 1.287 
5251 1.286 
5292 1.286 
Zo 1.285 
Bos 1.285 
587 5 izes 
aor 1.285 
5402 1.286 
5406 1.286 
5405 1.287 
5398 1.288 
B05 1.290 
5308 1.293 
5230 1.296 
5136 1.300 
SOsH 1.304 
4923 1.308 
4809 i. siz 
4694 3 
4581 1.32 
4467 


1.324 


ye All quantities in standard British units 
120 psi; T = 
62079 h’= 80,000 Et 
4) Exhaust data for frozen flow expansion 


2) P92 
>) eelta 


3200 R 


aD 


HTPB 


s2oo 
~ 304 
- 308 
soul 
vols 
woo 
ones 
~318 
—sl9 
oad 
~322 
2323 
~324 
~325 
32/7 
~328 
oe So 
~330 
SEY 
o332 
2355 
0334 
2335 
~ 336 


TABLE 2 
PEPCODE DATA FOR EQUILIBRIUM ADIABATIC COMBUSTION OF HTPB 


® is Y R Co Tath Ye 

Om ~0367 1.265 53. £02 326 4812 1.282 
0.6 .0440 ie oul 53.194 oor: 5041 L276 
Oa ~0514 M257 53.364 ~ 305 5229 1.272 
aya) 0567 Loo 53 612 go 5376 1.269 
0.85 .0624 tT. 2Do 53.7 66 ~341 5437 1.268 
0.9 ~0661 1.254 53 oa: 2342 5489 1.268 
0.95 .0697 1.254 54.137 2344 55 oL eeu 
iG .0734 1.254 54.353 ~345 5565 acior 
1.1 ~0807 1.254 54.848 -348 5606 1.267 
1.2 -0881 i250 55.404 350 5614 1.268 
ies .0954 1.256 56. Oo7 ~ Cae 5586 1.270 
1.35e 0991 ee 5 / 56,496 355 5559 1 eval 
1.4 -1028 1.258 56.748 . 356 5525 E2738 
1.45 .1064 P2509 S7eli9 soo 5484 1.274 
iS -1101 1.260) 57.505 - 358 5437 1.276 
se ~1174 1.263 58.293 “Joo S352 1.280 
ad ~1248 1.266 59.096 ~361 5216 1.284 
Eas ease: 1.269 59 . 93a: . 3618 5094 1.288 
1.9 1395 Le27e 60.749 ~ 365 4971 L.292 
2.0 - 1468 eS 61.568 366 4848 1.296 
2.1 ~1541 1.279 62.403 - 368 4725 1.300 
262 - 1615 1.282 63.191 ~ 369 4606 1.304 
Ze - 1688 1.286 63.994 » 380 4487 1.308 
Remarks: 


1) All quantities in standard British units 
620 psi; T = 3200 R 
6.0; h = 80,000 ft; isentropic inlet 
4) Exhaust data for frozen flow expansion 


2) P5 _ 
3) Mo = 


56 


= 3m 
.316 
Sou 
~325 
< 327 
.328 
Ashe 
.332 
» 3314 
a 7) 
~339 
~340 
~340 
341 
.342 
~343 
344 
~345 
346 
~347 
.348 
~349 
349 


TABLE 3 
PEPCODE DATA FOR EQUILIBRIUM ADIABATIC COMBUSTION OF HTPB 


$ f Y R Cp Tath Ye Coe 
0.6 0440 1.267 52.876 ei Pa: 3946 275 298 
e -0514 ire Ou 52.913 329 4218 1.288 ~ 304 
0.8 -0587 1.257 533035 2334 4436 1.282 a0 
0.9 -0661 1.254 53.295 530 4595 eas 2 314 
0 .0734 Lie oS 35601 3 341 4694 eor d sole 
ol 0807 1.254 54.183 2344 4735 1.278 ~ 320 
2 -O88l a2 52 54.826 0347 4718 a 2oo «322 
3 0954 1.258 55.570 ~ 348 4648 1.284 7o25 
1.4 » L028 Teor 56.389 23910 4546 1.289 0324 
5 lO 1.266 57 357 ao 4415 1.294 .324 
6 1174 Le2i69 58.098 soo e 4302 2299 2324 
dae? 1248 liZ7 3 58.956 wS3 4178 1.304 eo 
8 092 T e277 59.806 1354 4052 1.309 woz 
a) . 1895 ee 60.640 Do 3932 m 314 »326 
0 - 1468 1 #285 61.498 ~ 356 3810 1.319 «Sew 
J. ~1541 dee 289 62.555 oe 3690 1.324 woe] 
2.2 « EHD 1.294 63.144 +306 3576 oe -328 
Remarks: 


1) All quantities in standard British units 
2) Po = 50 psi; T = 1640 R 

3) Mo = 4.0; h = 80,000 ft 

4) Exhaust data for frozen flow expansion 
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APPENDIX B - RAMJET 
PROGRAM RAMJET 


Pete ete Hts eke Fee HK Fe KCL CLE BEA CT TELNAES TEAR ERIC REET TG 


ae a oe oe ee ee ee ee ee ee ee re a er oe ee ee a 


PROGRAM USED TO DETERMINE THE EXIT CONDITIONS OF A RAMJET BURNER ** 


CHAMBER. A LEGEND FOR THE VARIABLES INVOLVED IS OFFERED BELOW. 


STATION 0: FREESTREAM CONDITIONS 


STATION 2: INPUT CONDITIONS FOR THE COMBUSTION CHAMBER 


STATION 4: EXIT CONDITIONS OF COMBUSTION CHAMBER 
STATION E: EXHAUST CONDITIONS(POST NOZZLE) 


M: MACH NUMBER 

T: STATIC TEMPERATURE (DEG R) 

TT: TOTAL TEMP (DEG R) 

P: STATIC PRESSURE (LBF/IN2) 

PT: TOTAL PRESSURE (LBF/IN2) 

GAMMA: RATIO OF SPECIFIC HEATS 

R: GAS CONSTANT (BTU/LBM-DEG R) 
CP:SPECIFIC HEAT (BTU/LBM-DEG R) 

F: FUEL TO AIR RATIO 

PHI: EQUIVALENCE RATIO 

MDOT: MASS FLOW RATE (LBM/SEC) 

U: FLUID VELOCITY (FT/SEC) 

GC: PROPORTIONALITY CONSTANT (FT-LBF/LBM-SEC2) 
J: CONVERSION FACTOR = 778 BTU/FT-LBF 
DELTA: CONVERGENCE LIMIT 

D: PORT DIAMETER (IN) 

STEP: MACH NO. INCREMENT 

THRUST: UNINSTALLED THRUST (LBF) 

STHRST: SPECIFIC THRUST (LBF-SEC/LBM) 

SFC: SPECIFIC FUEL CONSUMPTION (LBM/LBF-SEC) 
ETA: THERMAL EFFICIENCY (TEMPERATURE RISE) 
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weve 
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Jove 


teve 
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DOUBLE PRECISION PO,P2,P4,PTO,PT2,PT4,T0,T2,T4,TTO,TT2,TT4,M2,M4 
DOUBLE PRECISION CP2,CP4,GAMMA2,GAMMA4G,R2,R4,F,H,MDOTO,GC,CD,U2 
DOUBLE PRECISION U4,TEST,DELTA,MDOT4,STEP,DELTT,MO,UE,U0, THRUST 
DOUBLE PRECISION STHRST,D2,D4,PI,GAMMAE,ME,TE,PE,PTE,TTE,ETA 
DOUBLE PRECISION TT4TH,J,CPE,PHI,SFC 

CHARACTER*1 ANS, ANSWER 

PRINT *,’ALL UNITS ARE STANDARD ENGLISH UNITS’ 

INITIAL DATA 

DATA MO,PO0,TO/6,.404,390/ 

DATA CD,GC,DELTA,STEP/1.5,32.174, .0001, .00001/ 

DATA M2,GAMMAO,RO,CPO,D2/.2000,1.40,53.3,.24,1.674/ 

DATA R2,CP2,GAMMA2/53.3,.28,1.33/ 

DATA J,PT2,MDOTO/778,146.72,1.00D0/ 

PE = PO 

PI=ASIN(1.D0)*2.D0 

A2= (PI/4.D0)*D2*D2 

PRINT *,’A2=',A2 

CONSTANT AREA COMBUSTOR 

A4=A2 

INITIAL FLOW COMPUTATIONS 
PTO=PO*(1.+MO*MO*(GAMMAO-1.)>*.5)%*(GAMMAO/(GAMMAO-1.)) 
TTO=TO*(1.+MO*MO*(GAMMAO-1.)*.5) 

TT2=TTO 

T2=TT2/(1.+M2*M2*(GAMMA2-1.)*.5) 
P2=PT2/((1.+M2*M2*(GAMMA2-1.)*.5)#*(GAMMA2/(GAMMA2-1.))) 
PRINT 5,’P2=',P2 

FORMAT(1X,A3,F9.4) 

PRINT 6,’T2=',T2 

FORMAT(1X,A3,F9.4) 

INITIAL GUESS ON EXIT MACH 

PRINT *, ‘INPUT ESTIMATE OF M4’ 

READ *, M4 

INPUT PEPCODE DATA 
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PRINT *,’INPUT TT4TH(DEG R)’ 
READ *, TT4TH 
PRINT *,’INPUT GAS PROPERTIES (F,R4,CP4,GAMMA4, GAMMAE) ’ 
READ *, F,R4,CP4,GAMMA4 ,GAMMAE 
PHI = F/0.0734D0 
DETERMINE ETA FOR HTPB 
IF (PHI .LE. 1) THEN 
ETA = 1.D0 - .35D0*PHI**.92D0 
ELSE 
ETA = .3D0 + .35D0*PHI**.92D0 
END IF 
SOLVE FOR M4 
TT4 = ETA*(TT4TH - TT2) + TT2 
IF (M4 .GE. .999) THEN 
PRINT *,’FLOW IS CHOKED FOR THIS VALUE OF F’ 
GOTO 50 
ENDIF 
CPE = R4*GAMMAE/((GAMMAE-1.D0)*J) 
T4=TT4/(1.+M4*M4* (GAMMA4-1.)*.5) 
U4=M4*SQRT (GC*GAMMA4*R4*T4) 
U2=M2*SQRT (GC*GAMMA2*R2*T2 ) 
MOMENTUM EQUATION FOR P4 
P4=P2+(MDOTO*(U2-(1.+F)*U4) )/(GC*A2) -CD*P2*U2*U2/(2.*R2*T2*GC) 
PT4=P4*(1.+M4*M4* (GAMMA4-1.)*.5)%*(GAMMA4/(GAMMA4-1.)) 
MDOT 4=MDOTO*(1.+F) 
TEST = MDOT4 - P4*A4*U4/(R4*T4) 
CHECK GUESS ON M4 
IF (ABS(TEST) .GE. DELTA) THEN 
M4=M4+STEP 
GOTO 10 
END IF 
SOLVE FOR EXIT CONDITIONS 
ME=SQRT((2.DO/(GAMMAE - 1.D0)*((PT4/PE)**((GAMMAE-1.D0)/GAMMAE) 
- 1.D0))) 
UE=SQRT(2.*GC*J*CPE*TT4* (1. -(PE/PT4)%**( (GAMMAE-1.) /GAMMAE) ) ) 
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20 


25 
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35 


40 


45 


60 


65 


70 


13> 


80 


85 


90 


TE=TT4/(1.DO+ME*ME* (GAMMAE-1.D0)*. 5D0) 
U0=MO*SQRT(GC*GAMMAO*R2*TO) 
THRUST=(MDOT4*UE -MDOTO*UO) /GC 
STHRST=THRUST/MDOTO 

SFC = F/STHRST 

PRINT 15,'MO’,’M2’, M4’, 'ME’ 
FORMAT(/1X,4A9) 

PRINT 20,MO,M2,M4,ME 

FORMAT(1X, 4F9.3/) 

PRINT 25,’T0’,'T2’,'T4", "TE! 
FORMAT(1X, 4A9) 

PRINT 30,T0,T2,T4,TE 
FORMAT(1X,4F9.1/) 

PRINT 35,/TTO’, TT2’, 'TT4’, *TTE’ 
FORMAT(1X, 4A9) 

PRINT 40,TTO,TT2,TT4,TT4 
FORMAT(1X,4F9.1/) 

PRINT 45,'PO0’,'P2’,’P4’, 'PE’ 
FORMAT (1X, 4A9) 

PRINT 60,P0,P2,P4,PE 
FORMAT(1X,4F9.3/) _ 

PRINT 65,’PTO’,’PT2’,’PT4’, PTE! 
FORMAT(1X, 449) 

PRINT 70,PTO,PT2,PT4,PT4 
FORMAT(1X, 4F9.3/) 

PRINT 75,'U0’,'U2', U4", "UE" 
FORMAT(1X, 4A9) 

PRINT 80,U0,U2,U4,UE 
FORMAT (1X, 4F9.3/) 

PRINT 85,’F’,’MDOTO’, THRUST’ , "SPTHRUST’ , ’SFC’ 
FORMAT(1X, 5A9) 

PRINT 90,F,MDOTO, THRUST , STHRST, SFC 
FORMAT(1X,4F9.4,F10.6/) 

WRITE (*,*) "PERFORM ANOTHER RUN THIS FLIGHT CONDITION (Y/N)?’ 
READ (*,’(A)’) ANS 
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IF (ANS .EQ. ’N’) GOTO 110 
GOTO 50 
110 WRITE (*,*) 'START NEW PROBLEM (Y/N)?’ 
READ (*,’(A)’) ANSWER 
IF (ANSWER .EQ. ’N’) STOP 
GOTO 100 
END 


SAMPLE OUTPUT: 


run raml fortran 
VS FORTRAN COMPILER ENTERED. 14:06:42 


*xKRAMJET*#* END OF COMPILATION 1] %ieieinik 
VS FORTRAN COMPILER EXITED. 14:06:45 


EXECUTION BEGINS... 

ALL UNITS ARE STANDARD ENGLISH UNITS 
A2= 2.20090103 

P2= 142.8812 

T2=3177.0289 

INPUT ESTIMATE OF M4 

? 


eo 
INPUT TT4THCDEG R) 
? 


5553 
aU GAS PROPERTIES (F,R4,CP4,GAMMA4 , GAMMAE) 
-0734,54.848, .344,1.258,1.285 


MO M2 M4 ME 
6.000 0.200 0.293 4.309 
TO T2 T4 TE 
390.0 3177.0 4548.4 1261.3 
TTO eo TT4 TTE 
3198.0 3198.0 4598.7 4598.7 
PO P2 P4 PE 
0.404 142.881 130.670 0.404 


PTO Peo PT4 PTE 
637.867 146.720 137.8795 1377s 


UO UZ U4 UE 
5805.821 538.372 931.7082 728722 


F MDOTO THRUST SPTHRUST SFC 
0.0734 1.0000 62.6804 62.6804 0.001171 


PERFORM ANOTHER RUN THIS FLIGHT CONDITION (Y/N)? 
n 

START NEW PROBLEM (Y/N)? 

n 

R; T=0.97/1.54 14:07:18 


62 


APPENDIX C - MIXER 
PROGRAM MIXER 


CHARACTER ANS 

REAL ETA,M4,T4TH,TT4TH,TT4,1T4,PT4,P4,A4,A3A,A3B, GAMMA4 ,R4,CP4 
REAL TT3A,TT3B, TT3BAR,U3A,U3B,U4 ,MDOT3A,MDOT3B,MDOT4 , GC 

REAL FNT,FNTT,FNU, FNA, FNGAM, FNP, P3A,P3B,V4,DELPA, DIFF 


C-—---INPOT ALL THE PERTINENT DATA 
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PRINT *, ‘INPUT R4,GAMMA4 ,T4TH!' 

READ *,R4,GAMMA4 , TATH 

PRINT *, ‘INPUT MDOT3A,MDOT3B,TT3A,TT3B' 
READ *, MDOT3A,MDOT3B,TT3A,TT3B 

PRINT *, ‘INPUT P3A,P3B,A3A,A3B! 

READ *,P3A,P3B,A3A,A3B 

PRINT *, ‘INPUT U3A,U3B' 

READ *,U3A,U3B 

PRINT *, ‘INPUT ETA,M4! 

READ *, ETA,M4 


GC = 32.174 

TT4TH = FNTT(T4TH, GAMMA ,M4) 

MDOT4 = MDOT3A+MDOT3B 

TT3BAR =(MDOT3A*TT3A + MDOT3B*IT3B) /MDOT4 
TT4 = (TT4TH - TT3BAR)*ETA + TT3BAR 

T4 = FNT(TT4,GAMMA4 M4) 

U4 = FNU(GAMMAG,R4,T4,M4) 


Ne 


A4 = A3A + A3B 

PT4 = (MDOT4*SORT(TT4) *FNA(GAMMA4 ,M4) ) / (A4*FNGAM(GAMMA4 , R4) ) 
P4 = FNP(M4,PT4,GAMMA4) 

DELPA = P3A*A3A + P3B*A3B - P4*A4 

V4 = (.9*GC*DELPA + MDOT3A*U3A + MDOT3B*U3B) /MDOT4 


DIFF = (U4 = V4) 


PRINT 15,TT3BAR 
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15 FORMAT(1X, ’TT3BAR =’,F9.3) 
PRINT 20,TT4TH 

20 FORMAT(1X,’TT4TH =’,F9.3) 
PRINT 25,TT4 


25 FORMAT (1X, ’TT4 =',F9.3) 
PRINT 30,T4 

30 FORMAT (1X, ’T4 =',F9.3) 
PRINT 35,04 

35 FORMAT (1X, ’U4 =' ,F9.3) 
PRINT 40,PT4 

40 FORMAT (1X, ’PT4 ae Se) 
PRINT 45,P4 

45 FORMAT (1X, ’P4 af F973) 
PRINT 50,V4 

50 FORMAT (1X, ’V4 =’ F9.3) 


PRINT 55,DIFF 
55 FORMAT(1X,’U4-V4 =" , F9.3) 

IF (ABS(DIFF) .LT. .01) STOP 

PRINT *, CONTINUE? (Y/N)’ 

READ (*,’(A)’) ANS 

IF (ANS .EQ. ’N’) STOP 

GOTO 10 

END 
ee PNT Tiececiciiccideciiciciciciicicicicicieiciciciccicicicicicicicicriricicseecsrsriecierkiciciesesrsrreieck 
ek ‘THIS FUNCTION DETERMINES STAGNATION TEMPERATURE (DEG R) GIVEN ex 
4k ~~ STATIC TEMPERATURE, MACH AND GAMMA. dele 
kkk RKC KR ;_e,C—_—RK—_—M—_ «M}$pM9 } PKR} Keiko aad am aR RgI aR ERIK aa: 

FUNCTION FNTT(T,GAM,M) 

REAL T,GAM,M 

FNTT= T*(1.+.5*M*M*(GAM - 1.)) 

END | 
toeiek BN Tice iceicdccicicciciciceicicicericiciecicicecieeececcscceiciersescsicsestsesessses 
#cek THIS FUNCTION DETERMINES STATIC TEMPERATURE GIVEN MACH, GAMMA *int 
#ckk AND STAGNATION TEMPERATURE. etese 
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KKEKKERE KKK tice kei“ cw clloclclelcloclecdelcteteiok 
FUNCTION FNT(TT,GAM,M) 
REAL TT, GAM, M 
FNT = TT/(1.+.5*M*M*(GAM - 1.)) 
END 
FFI CIS BN FCRIII CRI CRICRICKIERICIRIORIOREIRICRIRRERICRERRIARERIRIREEREREREREREERICR SG 
wie THIS FUNCTION DETERMINES GAS VELOCITY GIVEN GAMMA, GAS CONSTANT 2% 
wie STATIC TEMPERATURE AND MACH. devote 





FEV LCT L ECLECTIC. 


FUNCTION FNUC(GAM,R,T,M) 
REAL GAM,R,T,M,GC 


GC = 32.174 
FNU = M*SQRT(GC*R*GAMAT ) 
END 


KR KEN AK RK KK KE RII RAR RAE REI RARER AERA ERERERERERICKE 


Wick THIS FUNCTION DETERMINES THE RATIO A/A* GIVEN MACH AND GAMMA. *** 


RHF KH KERR EE CK KK RRR ASSERT AAAI IRAE RRR RR eee 


FUNCTION FNA(GAM,M) 


REAL GAM,M 

FNA = (1.0/M)*((2.0/(GAM+1.0))*(1.0+.5*M*M*(GAM - 1.0)))%*( (GAM 
+ +1.0)/(2.0*(GAM - 1.0))) 

END 


KKKKKENG AMOR ek nidcdncdccccicrciiidicdiccecdcdecdercicticcclccneciidecvdeckink 


wee THIS FUNCTION IS THE FUNCTION OF GAMMA AND R IN THE CONTINUITY **¥% 
wick EQUATION. ered 
Jenoininididinidinicindicinidicidvicieicicicicincicicicicicicicicicieicicicicicicicicicicicieicicieieleicieieeteieieteteiereseketetese 
FUNCTION FNGAM(GAM, R) 
REAL GAM,R,GC 


GC = 32.174 
FNGAM = SQRT((GAM*GC/R)*(2.0/(GAM+1.0))**( (GAM + 1.)/(GAM - 1.))) 
END 


Ye ve deve de He BN P ve deve Hote He ve Fe Fe KK LIA ALAA ALA LR AA RACAL ADA RTA AAAS SEDC LLNS REINER 
wick THIS FUNCTION DETERMINES THE STATIC PRESSURE GIVEN GAMMA, MACH ** 
IK AND STAGNATION PRESSURE. wee 


Hk Fete dete tek eke Ke SK Le AAA A AAA AAAI ATA ARAN AIA AKAIARII A AAR ARAMA HINA IRR ESERIES 
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FUNCTION FNP(M, PT, GAM) 

REAL M,GAM, PT 

FNP = PT/(1.0 + .5*M*M*(GAM - 1.0))**(GAM/(GAM - 1.)) 
END 


SAMPLE OUTPUT: 


run mixer fortran 
VS FORTRAN COMPILER ENTERED. 14:12:36 


*kMIXER** END OF COMPILATION 1 Aievederrs 
*kENTT*¥* END OF COMPILATION 2 *inivivk 
*kENT** END OF COMPILATION 3 *#itiekrex 
*xXENU** END OF COMPILATION 4 %itierei% 
**XFNA** END OF COMPILATION 5 inieininink 
**FNGAM** END OF COMPILATION 6 *tieekk 
*XENP** END OF COMPILATION 7 *#ieiekkk 
*KENPT** END OF COMPILATION 8 *ktierek 
VS FORTRAN COMPILER EXITED. 14:12:54 
EXECUTION BEGINS... 

ee R4 , GAMMA4 , TT4TH 


54.415,1.258 ,5292 
INPUT CP3A,CP3B,CP4 


t 
Z4EBG e502, 24k 
ele MDOT3A ,MDOT3B,TT3A,TT3B 


-9,.9661, 3198 , 4849.5 
INPUT P3A,P3B,A3A,A3B 


? 
16.45 0Gr4 ay 1 Os 
INPUT U3A,U3B 


? 
4666.61,5594.2 - 


ae M4 

1.423 

TT3BAR = 4193.020 
TT4TH = 5292.000 
TT4 = 5182.102 
T4 = 4108.812 
U4 = 4280.703 
PT4 = 106.573 
P4 = S497 2 
V4 = 4280.520 
U4-V4 = 0.184 
MDOT4* = 1.0660 
MDOT4 = 1. OGi62 


CONTINUE? (Y/N) 
n 
R; T=0.58/1.46 14:13:56 
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APPENDIX D - SFRJ Properties 
TABLE 4 


PROPERTIES THROUGH SFRJ 


® S20 ~10 1.00 30 io 0 
NAT .82 .75 .65 a75 .90 
Mo iy. Oo 1.00 1.00 fio oOo 
mf .04 .05 SON 5 1G a3 
T9 3L986:00 3198 300 31918. 00 3198.00 3198.00 
T4 4459.70 4611.40 4598.70 4838.00 4849.50 
Ma n27 28 n29 .32 .34 
Cy, 4416.80 4564.10 4548.40 4776.40 A7 73790 
Po 637.37 637.37 637.07 687.07 63/557 
Po 146.73 146.73 146.73 146.73 146.73 
Pig 138.79 £35732 37789 136.93 35,70 
Po .40 .40 Ta .40 .40 
P2 142.89 142.89 142.89 142.89 142.89 
Pa m2. 56 Sass 130.68 128.65 zor. O7 
Uo 5a05. 80 Ba Oo 0 5S05—780 5805.80 58 0ER50 
U5 SSies 7 Scone 7 5aie > 7 536.37 536437 
Ug Ssioe7 7 889.56 930598 1043.38 1170. oS 
Ue OSs 3o.1 7201.00 7287.61 7566.49 7727.80 
F 45.59 54.86 62.68 77 <XO 91.47 
S - 8.0492e-4 9.3650e-4 1.1710e-3 1.2366e-3 1.4444e=-3 
Remarks: 

1) Mo = 6.0; h = 80,000 ft 

2) See = in( o) 

3) All quantites in standard British units 
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TABLE 5 


PROPERTIES THROUGH SFRJ 


® 220 70 1.00 reo 1.50 
NAT 82 a7 5 .65 .75 .81 
Mo IO 1.00 1.00 i200 1.00 
mM .04 .05 SOY 0 al 
T9 3198.00 3198.00 3198.00 3198.00 3198.00 
T4 4513.40 4716.80 4736.50 4978.20 5007.70 
My, ei 029 ~ 30 o2 ao 
ty, 4469.80 4668.00 4684.10 4913.70 4936.80 
Po 637.87 637.87 637.87 637.87 637.87 
Po 637.87 637.87 "63/.87 637.87 637.87 
P4 603.09 600.82 598.77 594.47 DyZ. LO 
Po .40 40 40 40 40 
Po 621.18 621.18 621.18 621.18 621.18 
P4 57.75 571.00 566.72 557.64 552.74 
ug 5805.80 5805.80 5805.80 5805.80 5805.80 
ug 53605 / 538.37 5338.37 536n 3 / DOO. 37 
ug 843.85 905.64 952.10 1067.93 1 25607 2 
Ue 7487.60 7736.89 7857.18 OLS ono Ss 8244.19 
iF 60.81 72.38 81.68 97.25 104.00 
S 6.0349e-4 7.099le=-4 8.9859e-4 9.8122e-4 1.0587e=-3 
Remarks: 

1) Mp = 6.0; h = 80,000 ft 

2) nat = En( & ) 

3) All quantites in standard British units 

4) Isentropic inlet 
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® ~ 70 
NAT wn 
Mo 1.00 
M+ 205 
T> 1638.00 
T, 3567.40 
M4 37 
es 3504.30 
Re 61.34 
P5 31915 
Py 29.26 
Po ~404 
P2 31.09 
P4 26.c5 
Uo 3870.55 
us 390.90 
ug BOs. 35 
Ue 5785.54 
68.76 
S 

Remarks: 


~90 

- 68 
1.00 
«07 
1638.00 
3655.10 
~o2 
8586.50 
61.34 
31.95 
2014 
~404 
31.09 
26.01 
3870.55 
390.910 
1078.20 
5904.23 
J. 33 


TABLE 6 
PROPERTIES THROUGH 


1.10 

- 68 
1.00 
08 
1638.00 
3749.80 
241 
3672.10 
61.34 
31.95 
BO «92 
~404 
S102 
26.07 
3870.55 
390.90 
1656 «89 
6030.25 
82.26 


Mo = 4.0; h = 80,000 ft. 
nat = fn( © ) 


.All quantites in standard British units 
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SFRJ 
E30 
so 
1.00 
~ LO 
1638.00 
Bool.ou 
044 
o789 .20 
61.34 
31.95 
28.63 
~404 
Oe 
25.45 
3670.55 
390.90 
fa/1.23 
6188.78 
90.41 


1.80 
~90 
1.00 
oS 
1638.00 
3813.00 
2 48 
3696.50 
61.34 
Sino.) 
28.18 
~404 
SA. 62 
24.42 
S670.55 
390.90 
c437.11 
6272.16 
100.40 


7.4725e-4 8.7692e-4 9.8154e-4 1.0554e-3 1.3159e=-3 


APPENDIX E - DMRJ Properties 
TABLE 7 
PROPERTIES THROUGH DMRJ (50/50) 


Poppe 1.20 1.50 1.80 ZO 2520 
O70T -60 -75 -90 1.10 1.25 
M3a 50 .50 .50 -50 50 
M37 .54 56 57 .58 .58 
Tan 3198.00 3198.00 3198.00 3198.00 3198.00 
T3B 4774.50 4903.40 4849.50 4612.80 4534.30 
Tg 4856.42 5056.25 5182.10 5268.11 5265 . 50 
M3p 2.06 2.05 2.05 2.04 2.04 
M4 1.54 1.49 1.42 L 30 1.27 
t3p 2976.70 3034.30 2922.50 2766.30 2706.50 
C4 3701.93 3927.17 4108.81 4312.30 4359.34 
Po 637.87 637.87 637.87 637.87 637.87 
Poa 281.96 281.96 281.96 281.96 281.96 
Pop 146.73 146.73 146.73 146.73 146.73 
Daa 281.96 281.96 281.96 281.96 281.96 
P3p lic. OL 136.36 35.72 135. 17 5.09 
Py Is . 06 109.98 106.57 102.67 101.72 
Po - 40 - 40 -40 40 - 40 
P3a 16.40 16.40 16.40 16.40 16.40 
P3B 16.40 16.40 16.40 16.40 16.40 
P4 30.70 32.20 34.357 39.04 40.58 
U3z 4666.61 4666.61 4666.61 4666.61 4666.61 
U3B 5391.56 5539.50 5594.20 5568.70 5548.60 
Ug 4365.81 4630.70 4280.70 4052.73 3971.72 
Ue 7295.48 7488.00 7620.32 7730.24 7748.30 
Ug 5805.80 5805.80 5805.80 5805.80 5805.80 
F 56.28 65.11 72.05 79.20 80.70 
S 7.8183e-4 8.4626e-4 9.1738e-4 1.0194e%3 1.0458e-3 
Remarks: 


1) Mg = 6.0; h ©980,000 ft 


\O 
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TABLE 8 
PROPERTIES THROUGH DMRJ (50/50) 


@srprg ee 0 ie 0 1.80 2220 2.50 
Sror 60 275 .90 ie © iene 
M30 .50 .50 .50 .50 0 
N3p .54 .56 oi 58 .58 
T3a S95.0 0 3198.00 3198.00 3198.00 3198.00 
T3B 4774.50 4903.40 4849.50 4612.80 4534.30 
Ty, 4611.70 4774.99 4852.41 4859.43 4848.02 
M3, 2.06 2 305 an05 2.04 2.04 
M4 1.69 1.66 lez io) i 53 
t35 2976.70 3034.30 Zo2e5 50 2766.30 2706.50 
4 Sie 52. 56 3523.64 SiGzowos 3708.57 3721.38 
Po 637.87 637.87 637 797 637.87 637.87 
Poa 281.96 281.96 281.96 281.96 281.96 
Pop 146.73 146.73 146.73 146.73 146.73 
P3a3a 281.96 281.96 281.96 251996 281.96 
P3B lest cO 136.36 1SiS7 2 igo 17 i135. 09 
Py 1235s ee 120.22 dees 7 ies. OF 112.69 
Po 40 ~40 -40 -40 40 
P3a 16.40 16.40 16.40 16.40 16.40 
P3B 16.40 16.40 16.40 16.40 16.40 
P4 26.63 27.44 28.43 30.44 Sal Os 
U3a 4666.61 4666.61 4666.61 4666.61 4666.61 
U3B Seo S16 553 or 50 5594.20 5568.70 5548.60 
U4 4559.37 4590.77 45710. 16 4470.29 4433.61 
Ue 7s 5.60 7305.35 7405.50 7465.70 7470.70 
Ug 5805.80 5805.80 5805.80 5805.80 5805.80 

bale 09 Some 2 64.93 70.3 3 71.34 
S 8.6123e-4 9.3203e-4 1.0179e-3 1.148le-3 1.1830e-3 
Remarks: 


1) Mo = 6.0; h = 80,000 ft 


ad 


TABLE 9 
PROPERTIES THROUGH DMRJ (50/50) 


Popry ES 2x0 1.50 1.80 2.20 2.30 
dror 60 as -90 sO aL AE 
M3, a0 a0 a0 oo0 50 
3p 54 .56 7 .58 .58 
T3a 3198.00 3198.00 3198.00 3198.00 3198.00 
Ts, 4774.50 4903.40 4849.50 4612.80 4534.30 
T4 4530.13 4681.24 4742.51 | 4726.53 4709.53 
M35 2.06 2505 2-05 2.04 2.04 
M4 1.74 digeee? 1. 1.68 Linas 1.62 
t3p 2976.70 3034.30 2922.50 2766.30 2706.50 
Cy, 3239.07 3392.20 3473.37 3518.96 35922676 
eG 6s7o7 637 4157 627 .8/ 637.87 657.07 
Pon 281.96 281.96 281.96 281.96 281.96 
Pop 146.73 146.73 DEN FS 146.73 146.73 
Pan 281.96 281.96 281.96 2815296 281.96 
P3p 138.01 136.36 17 2 125ea7 135.09 
Py 127 22 1457 122.00 118.82 117.90 
Po 40 - 40 40 ~ 40 ~ 40 
P3a 16.40 16.40 16.40 16.40 16.40 
P3B 16.40 16.40 16.40 16.40 16.40 
P4 25.40 26.103 26s 2 28.19 28.62 
Uz, 4666.61 4666.61 4666.61 4666.61 4666.61 
U3B 539 17..56 5539.50 5594.20 5568.70 5548.60 
Ug 4616.78 4659.38 4654.91 4579.11 4550.24 
Ue 7082.10 7244.40 7333480 T3772 38 7375.00 
Uo 5805.80 5805.80 5805.80 5805.80 5805.80 

Aor 0 5 JZ 62.56 67.36 68.22 
S 8.9152e-4 9.6464e-4 1.0566e-3 1.1987e-3 1.2372e-3 
Remarks: 


i) WMG. = 96.0; Hijemeenooo aur 
2 ) NAT = 0.5 
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?SFRI 
fTOT 
M3A 
M3B 


P3A 
P3B 
P4 

UZA 
U3B 


Remarks: 


1) 
2) 


1.20 

- 60 

~50 

~54 
3198.00 
4774.50 
4326.19 
2.06 
1.88 
2976.70 
2957.48 
637.07 
281.96 
146.73 
265.90 
SS oo"u 
Boo 395 
-40 
16.40 
16.40 
Las 
4666.61 
So 7ie oS 
a7 SSO 
6948.42 
5805.80 
45.02 


TABLE 10 
PROPERTIES THROUGH DMRJ (50/50) 


350 
a1 5 
~50 
~56 

3198.00 
4903.40 
4446.86 
2.05 
1.86 
3034.30 
3068.96 
637257 
251.96 
146.73 
281.96 
Ea6eroo 
raise 
~40 
16.40 
16.40 
22.78 
4666.61 
5539.50 
Aelia 
NOOR 2 
5805.80 
52 oe 


my = 6.0; Hh ="so;000 ft 
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1.80 
90 

.50 

noi 
3198.00 
4849.50 
4467.76 
2.05 
1.86 
2922.50 
3094.25 
637.87 
281.96 
146.73 
281.96 
13S5e72 
ie26 
40 
16.40 
16.40 


2276S 


4666.61 
5594.20 
4842.54 
7153 228 
5805.80 

56.60 


2.20 
1.10 
~50 

-58 
3198.00 
4612.80 
4394.30 
2.04 
1.84 
2766.30 
3058.51 
637.07 
281.96 
146.73 
281.96 
1357, 
136.36 
~ 40 
16.40 
16.40 
2a ea 
4666.61 
5568.70 
4816.07 
T2554. 51 
5805.80 
S9e oT 


9.7743e~-4 1.0573e-3 1.1679e-3 1.3478e-3 1l. 


2 sa 
das 
50 

58 
3198.00 
4534.30 
Aaos 30 
2.04 
1.84 
2706.50 
3040.95 
So7ee7 
281.96 
146.73 
281.96 
135.09 
eee 
40 
16.40 
16.40 
eel 
Pecaeen 
5548.60 
4803.16 
7145.50 
5805.80 
60.38 
3977e-3 


TABLE 11 
PROPERTIES THROUGH DMRJ (80/20) 


Pcorri Leo eo LeSo 1.80 2.20 
707 88 1.04 1.20 1.44 1.76 
M3a ~20 w20 e20 © 20 ~20 
M3, 86 -88 89 2 ~93 
T3A 3198.00 3198.00 3198.00 3198.00 3198.00 
T3B 4693.40 4838.00 4903.40 4849.50 4612.80 
T4 5194.71 5293). 08! 5330.96 5271.16 5007 oS 
M33 2607 205 2205 2.05 2.04 
M4 L296 1.58 1.60 1.61 1,62 
t3p 2928.40 3010720 3034.30 2963.10 2753.06 
ty, 3958.40 4001.22 4009.98 3936.22 3714.24 
Po 637.87 637.87 637.87 637687 637.87 
Poa 281.96 281.96 281.96 281.96 281.96 
Pop 146.73 146.73 146273 146.73 146.73 
P3a 281.96 281.96 281.96 281.96 281.96 
P3B 137-56 136.92 136.35 13 Sam 13 506 
Py 103.60 104.79 1GS.60 106343 EO7 150 
Po -40 -40 -40 -40 -40 
P3a 16.40 16.40 16.40 16.40 16.40 
P3B 16.40 16.40 16.40 16.40 16.40 
P4 2) oS 26.78 26.35 26.07 26.04 
U3A 4666.61 4666.61 4666.61 4666.61 4666.61 
U3B 5323.47. 5450.17 5539.40 5594.10 5562.50 
U4 4590.69 4723.43 4812.61 4873.04 4859.37 
Ue 7607.70 7738.90 7823.28 7856.59 77 sOmilz 
ug 5805.80 5805.80 5805.80 5805.80 5805.80 
iigzeo 78.44 84.13 Soo 91.06 
S 9.0629e-4 9.7278e-4 1.0472e-3 1.1803e-3 1.4188e=-3 
2 
Remarks: 


1) Mop = 6.0; h = 80,000 ft 
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P3A 
P3B 
P4 
UZA 
U3B 
Ug 
Ue 


Uo 


Remarks: 


1) 
2) 


TABLE 
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PROPERTIES THROUGH DMRJ (80/20) 


LO 
-88 
~20 
- 86 

3198.00 
4693.40 
4952.03 

2.06 

1.69 

2928.40 
3614.14 
657.87 
281.96 
146.73 
201.96 
137..56 
112.21 
~40 
16.40 
16.40 
24.16 
4666.61 
5323.47 
4975.56 
7454.65 
5805.80 
66.22 


o-oo 
1.04 
~20 

- 88 
3198.00 
4838.00 
5060.26 
2.05 
eet 1 
3010.20 
367 ae 
637.87 
281.96 
146.73 
281.96 
136.92 
Poe LS 
~40 
16.40 
16.40 
23.73 
4666.61 
5450.17 
4894.66 
Togs OS 
SS Uomo 
T3-.00 


6m0> h = 780,000" E£t 
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a0 
1.20 
~20 

89 
34938 ..00 


4903.40 


DLOD. 70 
Z.05 
1.72 

3034.30 

3693.04 

637.87 
281.96 
146.73 
281.96 
LS 63S 
113.84 
~40 
16.40 
16.40 
23.46 

4666.61 

5539.40 

4976.88 

7682.08 

5805.80 

13S 


v2 7000e—4 1,.0373e=3 1.1102e-3 1. 


1.80 
1.44 
~20 

on 
3198.00 
4849.50 
5052-56 
2.05 
aNerseS 
2963.10 
3628.42 
637.87 
294.96 
146.73 
2c 
1357-1 
114.66 
~40 
16.40 
16.40 
23.27 
4666.61 
5594.10 
5035 155 
VIAT S90 
5805.80 
84.77 


2520 
1.76 
210 

.93 
3198.00 
4612.80 
4811.14 
DOr 
1.73 
2753.00 
3435.36 
637.87 
281.96 
146.73 
281.96 
135.16 
115.48 
.40 
ESA 
16.40 
23.36 
4666.61 
5562.50 
5011.57 
7606.55 
5805.80 
86.51 


2469e-3 1.4934e=-3 


®SFRI 
SToOT 
M3A 
M35 


P3A 
P3B 
P4 

UZA 
U3B 


Remarks: 


1) 
2) 


TABLE 
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PROPERTIES THROUGH DMRJ (80/20) 


1.10 
-88 
~20 
86 

3198.00 
4693.40 
4871.13 

2.06 

1.74 

2928.40 
3501.87 
637.87 
281.96 
146.73 
281.96 
L3W7 56 
115.68 
~ 40 
16.40 
16.40 
23.14 
4666.61 
5323.47 
4831.23 


7403.50 - 


5805.80 
64.52 


12430 
104 
20 

.88 
3198.00 
4838.00 
4982.67 
2.05 
1.76 
3010.20 
3565.87 
637.87 
281.96 
146.73 
281.96 
136.92 
116.47 
40 
16.40 
16.40 
De a 
4666.61 
5450.17 
4946.70 
7544.30 
5805.80 
71.93 


1.50 
1.20 
~20 

89 
3198.00 
4903.40 
5030.63 
2.05 
1.77 
3034.30 
3588.53 
637.87 
281.96 
146.73 
281.96 
136.55 
Le G 
~40 
16.40 
16.40 
22:57 
4666.61 
5539.40 
5028.47 
7635.10 
5805.80 
77.76 


1.0012e-3 1.0608e-3 1.1329e-3 


Mo = 


6.0; h = 80,000 ft 


76 


1.80 
1.44 
20 

91 
3198.00 
4849.50 
4979.82 
2.05 
yr 
2963.10 
3527.28 
637.87 
281.96 
146.73 
281.96 
i Sree 
117.93 
40 
16.40 
16.40 
22.40 
4666.61 


5594.10 > 


5083.15 
7671.05 
5805.80 

83.18 


2.20 
176 
20 

93 
3198.00 
4612.80 
4745.50 
2.04 
1:77 
2753.00 
3344.93 
637.87 
281.96 
146.73 
281.96 
Tes elt 
118.50 
40 
16.40 
16.40 
22.54 
4666.61 
5562.50 
5056.53 
7562.83 


5805.80 


84.98 


1.2708e-3 1.5204e-3 


®SFRI 
°ToOT 
M3A 
M3B 


 ~+BOA 
P3B 
P4 

U3ZA 
U3B 


Remarks: 


1) 
2) 


1.10 
-88 
~20 
- 86 

3198.00 
4693.40 
4668.90 

2.06 

1.86 

2928.40 
i223. Ou 
63: 71287 
281.96 
146.73 
2e1.96 
137.56 
126.18 
~ 40 
16.40 
16.40 
20.74 
4666.61 
5323.47 
4963.02 
7275.70 
5805.80 
60.29 


TABLE 
PROPERTIES THROUGH DMRJ (80/20) 


ie oO 
1.04 
~20 

-88 
3198.00 
4838.00 
4788.71 
2.05 
Poo? 
3010.20 
3300.05 
637.87 
Zoo oO 
146.73 
281.96 
136,92 
126.41 
~ 40 
16.40 
16.40 
20.58 
4666.61 
5450.17 
5070.58 
7422.45 
5805.80 
67.85 


Mp = 6.0; h = 80,000 ft 


en 


14 


1.50 
1.20 
~20 

~ 89 
3495.00 
4903.40 
4842.95 
2.05 
1.88 
3034.30 
3331.88 
637.87 
281.96 
146.73 
ZO i7 6 
136.35 
E2677 
~ 40 
16.40 
16.40 
20.47 
4666.61 
5539.40 
5147.29 


To? 00 


5805.80 
mae 


1.80 
1.44 
~20 

eel 
3198200 
4849.50 
4797.74 
2.95 
1.88 
2963.10 
32776 
637.87 
281.96 
146.73 
281.96 
Loom 
127.55 
~40 
16.40 
16.40 
20.34 
4666.61 
5594.10 
5200.83 
7554.90 
5805.80 
ra ee) 


iO miae—3 121246e-3.1.1943es3 1.3349e-3.1. 


2.20 
1e76 
~20 

~93 
3198.00 
4612.80 
4581.75 
2.04 
1.88 
2753.00 
3113.25 
On 760 
281.96 
146.73 
281.96 
LS52.26 
227562 
~ 40 
16.40 
16.40 
20.55 
4666.61 
5562.50 
5168.88 
7454.49 
5805.80 
Sis 
poeroe=s 


TABLE 15 
PROPERTIES THROUGH DMRJ (80/20) 


eee 1.20 1.50 1.80 2.00 2 120 
Smo 96 1.20 1.44 1.60 1.76 
Tete 20 20 20 20 126 
M35 87 89 91 Ro 93 
ites 1638.00 1638.00 1638.00 1638.00 1638.00 
Tas 3837.10 3882.50 3813.00 ~ a72em0cmmeazOmcO 
Ty 4540.29 4602.13 4394.83 4215.56 4027.05 
M35 2.26 2.25 2.25 225 2225 
My was a3 Wee 1.47 1.53 
en 2236.70 2222.40 2140.00 2065.10 1979.60 
ty 3718.25 3754.16 3481.94 3262.25 3051.32 
Po 61.34 61.34 61.34 61.34 61.34 
ae Ai 77 41.77 Ag Alea A ea7 
Pop 31.95 31.95 31.95 31.95 31.95 
ae 41.77 Acer AMedey 41.77 eb Fi 
Ee 28.76 28.42 28.18 28.05 27.95 
Py 15.94 16.05 16.61 17.09 17.56 
Po 40 40 40 40 40 
Baa 2.44 2.44 2.44 2.44 2.44 
Bae 2.44 2.44 2.44 2.44 2.44 
Py 5.93 5.89 5.43 5.10 4.85 
Berl 3316.50 3316.50 3316.50 3316.50 3316.50 
aoe 5080.44 5193.07 5222.49 5216.18 5189.83 
U4 3744.51 3837.09 4002.05 4100.34 4160.59 
ae 6297.70 6405.80 6359.22 6284.75 6205.44 
Up 3870.55 3870.55 3870.55 3870.55 3870.55 
F 89.24 96.34 98.24 97.97 97.49 
S 7.9002e-4 9.1448e-4 1.0759e-3 1.1983e-3 1.3253e-3 
Remarks: 


1) Mp = 4.0; h = 80,000.ft 


78 


TABLE 


16 


PROPERTIES THROUGH DMRJ (80/20) 


ScopRrg 20 1350 1.80 2200 2.20 
707 | ~96 20 1.44 1.60 L276 
M3. .20 .20 .20 .20 "eo 
3p 87 .89 .91 .92 =o 
T3a 1638.00 1638.00 1638.00 1638.00 L638. 00 
T3B 3837.10 3882.50 3813.00 3728.00 3620.60 
T4 AS97795 4254.53 4098.51 3956.26 3801.00 
M38 2.26 2.25 2.25 2.25 2.25 
M4 1236 1.57 — y. 62 1.66 70 
t3p 2236.70 2222.40 2140.00 2065.10 1979.60 
ty 3206.91 3236.18 3046.01 2884.58 2723.74 
Po 61.34 61.34 61.34 61.34 61.34 
Poa 41.77 41.77 41.77 41.77 41.77 
Pop 38, 95 Jn. 9S 3.95 31%. 95 31.95 
Pan An. 77 Any. 7 7 Anna 7 A 7 7 ane 7 
P33 26.76 28.42 28.18 28.05 27.95 
Py £7. 69 Peo 1 18.40 18.87 Loes3 
Po ~40 ~40 ~40 ~ 40 ~40 
P3A 2.44 2.44 2.44 2.44 2.44 
P3B 2.44 204 a4 2.44 2.44 
P4 4.66 4.63 4.42 4.25 4.12 
U3Zz 336.50 3oul6.50 3316.50 3316.50 3316.50 
U3B 5080.44 5193.07 5222.49 5216.18 5189.83 
Ug 4111.46 4204.95 4297.20 4347.46 4374770 
Ue 6110.30 6214.80 6194.14 6138.50 607587 
Uo 2.57.0). 55 3870.55 3870.55 38720555 S370. 55 

83.00 89.88 92.97 92.89 92.93 
S) 8.4937e-4 9.8020e-4 1.1419e-3 1.2639e-3 1.3903e-3 
Remarks: 


1) Mp = 4.0; h = 80,000 ft 


72 


TABLE 17 
PROPERTIES THROUGH DMRJ (75/25) 


Poprg iene 1236 1.50 1.80 220 
nor S25 .975 Tap 1.35 1.65 
M32 .25 75 25 25 .25 
3p .81 82 .83 85 .87 
T3. 3198.00 3198.00 3198.00 3198.00 3198.00 
T3B 4693.40 4838.00 4903.40 4849.50 4612.80 
1, 5144.86 5257.99 5315.18 5300.77 5101.19 
M3p 2.06 2.05 2.05 2205 2.04 
Ma 1.55 1.56 1.57 1.55 1.55 
t3p 2928.80 3010.60 3034.90 2963.70 2767.00 
ta 3930.61 4000.49 4034.09 4035.50 3864.15 
Po 637.87 637.87 637.87 637.87 637.87 
Pon 281.96 281.96 281.96 281.96 281.96 
Pop 146.74 146.74 146.74 146.74 146.74 
Pan 281.96 281.96 281.96 281.96 281.96 
P3p 137.48 136.83 136.25 135.58 135.02 
Py 104.58 105.17 105.50 104.63 104.93 
Po 40 40. 40 40 40 
P3a 16.40 16.40 16.40 16.40 16.40 
P3R 16.40 16.40 16.40 16.40 16.40 
P4 28.20 27.74 27 27.91 28.04 
U3za 4666.61 4666.61 4666.61 4666.61 4666.61 
U3R 5322.90 5449.47 5538.60 5593.17 5567.56 
U4 4536.71 4647.34 4722.08 4733.79 4721.40 
ue 7552.75 7696.37 7793.20 7849.98 7826.360 
Uo 5805.80 5805.80 5805.80 5805.80 5805.80 

68.52 75.89 81.78 Syd: 92.26 
S 8.8438e-4 9.4349e-4 1.010le-3 1.1298e-3 1.3126e-3 
Remarks: 


1) Mp = 6.0; h = 80,000 ft 


80 


®SFRI 
oT? 
M3A 
M3B 


P3A 
P3B 
P4 

U3A 


Remarks: 


1) 
2) 


TABLE 


18 


PROPERTIES THROUGH DMRJ (60/40) 


1.30 
~78 

~ 40 

~ 66 
3198.00 
4838.00 
DOs. 59 
2.05 
1.52 
3010.20 
3927. 36 
637567 
281.96 
146.74 
2m, 916 
136.92 
108.15 
~ 40 
16.40 
16.40 
30.12 
4666.61 
5450.16 
4458.84 
7518.44 
5805.80 
66.62 


Mo = 6.0; 


1.50 
~90 

~ 40 

~ 67 
3198.00 
4903.40 
5199.88 
220) 
1.51 
3034.30 
4024.09 
OStic7 
281.96 
146.74 
262.96 
136.35 
106.97 
~40 
16.40 
16.40 
30.65 
4666.61 
5539.41 
4480.46 
7634.66 
5805.80 
72.53 


80,000 ft 


oul 


1.80 
POS 
~40 

68 
5196". 00 
4849.50 
5281.00 
2.05 
1.46 
2963.10 
4414.48 
637.87 
281.96 
146.74 
281.96 
135.71 
104.78 
~40 
16.40 
16.40 
32.14 
4666.61 
5594.11 
4438.38 


7742.18 


5805.80 
Poaag 


Dee 0 
Tare 
.40 

70 
3198.00 
4612.80 
5271.89 
2.04 

1. 37 
2766.30 
4239.85 
637.87 
281.96 
146.74 
281.96 
135.16 
101.54 
.40 
16.40 
16.40 
35.21 
4666.61 
5568.60 
4274.51 
7800.58 
5805.80 
85.52 


2.30 
1.38 
~40 

~ 70 
S196, 00 
4534.30 
5249.18 
2.04 
35 
2706.50 
4240.61 
637.87 
281.96 
146.74 
250596 
135.06 
100.88 
~40 
16.40 
16.40 
22690 
4666.61 
5548.52 
4227.20 
7oto.9> 
5805.80 
87.09 


8.6010e-4 9.1137e-4 1.0004e-3 1.1343e-3 1.1632e-3 


PSFRI 
®TOT 
M3A 
M3B 


P3A 
P3B 
P4 

U3A 
U3B 


Remarks: 


1) 
2) 
3) 


TABLE 
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PROPERTIES THROUGH DMRJ (50/50) 


1.20 

- 60 

~ 50 

~54 
3198.00 
4923.00 
4959527 
2.44 
1.74 
2737.00 
IDIGh 25 
637.87 
637.87 
637.87 
637.87 
596.23 
343.50 
~ 40 
37.08 
37.08 
69.08 
4666.61 
6072.34 
4813.48 
7760.03 
5805.80 
71.38 


Isentropic inlet 


1.50 
75 
~50 
56 

3198.00 
5007.70 
5206.12 
2.44 
1.66 
2749.70 
3849.85 
637.87 
637.87 
637.87 
637.87 
592.67 
327.83 
~40 
37 . 08 
370s 
74.57 
4666.61 
6215.17 
4784.97 
8000.26 
5805.80 
81.91 


Mo = 6.07 h = 80,000 ft 


82 


1.80 
~90 

~ 50 

~57 
3198.00 
4906.30 
5364.06 
2.44 

es OW 
2642.50 
4087.70 
637.87 
637.87 
637.87 
637.87 
590.18 
312.65 
~ 40 
37.08 
37.08 
81.74 
4666.61 
6248.42 
4677.20 
8152). 93 
5805.80 


89.73. 
6.1786e-4 6.7271e-4 7.3662e-4 7.7816e-4 8.2716e-4 


2.00 
1.00 

~ 50 

~57 
3198.00 
4785.70 
5425.81 
2.44 
one 
2544.80 
4202.44 
637.107 
637.87 
637.87 
637.87 
588.97 
303.03 
~ 40 
37.08 
37.08 
85.83 
4666.61 
6234.66 
4596.10 
8236.10 
5805.80 
94.33 


2.20 

1. 20 
.50 

. 58 
3198.00 
4638.40 
5454.29 
2nd 
1.45 
2435.30 
4307.33 
637.87 
637.87 
637.87 
637.87 
588.06 
297.10 
- 40 
37.08 
277. O08 
92.62 
4666.61 
6198.93 
4470.48 
8276.85 
5805.80 
97.56 
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